Kertesz Line in the Three-Dimensional Compact U(l) Lattice Higgs Model 



in 
o 
o 



i 

CO 

c3 



I 

C 

o 
o 



> 

On 
OS 

in 
m 
O 
in 
o 

03 



i 

c 

o 
o 



X 
S3 



Sandro Wenzel, Elmar Bittner, Wolfhard Janke, Adriaan M.J. Schakel, and Arwed Schiller 
Institut fiir Theoretische Physik, Universitat Leipzig, 
Augustusplatz 10/11, D-04109 Leipzig, Germany 
(Dated: 2nd February 2008) 

The three-dimensional lattice Higgs model with compact U(l) gauge symmetry and unit charge is 
investigated by means of Monte Carlo simulations. The full model with fluctuating Higgs amplitude 
is simulated, and both energy as well as topological observables are measured. The data show a 
Higgs and a confined phase separated by a well-defined phase boundary, which is argued to be 
caused by proliferating vortices. For fixed gauge coupling, the phase boundary consists of a line 
of first-order phase transitions at small Higgs self-coupling, ending at a critical point. The phase 
boundary then continues as a Kertesz line across which thermodynamic quantities are nonsingular. 
Symmetry arguments are given to support these findings. 



Being one of the few well-understood theories exhibit- 
ing charge confinement, the three-dimensional (3D) pure 
compact U(l) gauge theory plays a central role in the 
study of deconfinement transitions . An intriguing ex- 
tension obtains by coupling a scalar matter field to this 
confining gauge theory 0. The resulting Higgs model 
and its extensions have recently attracted considerable 
attention also in the condensed matter community as ef- 
fective descriptions of quantum critical phenomena Q. 
In a seminal paper, Fradkin and Shenker j3] studied the 
phase diagram of the model in the London limit, where 
the Higgs field has a fixed amplitude. They concluded 
that in the case of a Higgs field carrying one unit charge 
q = 1, it is always possible to move from the Higgs re- 
gion into the confined region without encountering sin- 
gularities in local gauge- invariant observables. As for the 
liquid-vapor transition, this is commonly interpreted as 
implying that the two ground states do not constitute 
distinct phases. This is supported further by symmetry 
considerations In the absence of matter fields, the 
3D pure noncompact U(l) gauge theory is characterized 
by a global magnetic U(l) symmetry. When magnetic 
monopoles are introduced, which are pointlike instanton 
solutions of the compact U(l) gauge theory, this global 
symmetry becomes anomalous and only the discrete sub- 
group Z of integer numbers survives. When the compact 
gauge theory is subsequently coupled to a Higgs field car- 
rying charge q, the magnetic symmetry is further reduced 
to the finite cyclic subgroup Z g of q elements. For q = 2, 
this recovers the known result that the model undergoes 
a continuous phase transition belonging to the 3D Ising 
universality class |2J, |5(. For q = 1, this argument ex- 
cludes a continuous phase transition because the group 
Z±, consisting of only the unit element, cannot be spon- 
taneously broken. 

In this Letter, we argue that Monte Carlo data show 
a more refined picture with two distinct phases sepa- 
rated by a well-defined phase boundary. We consider 
the q = 1 model with fluctuating Higgs amplitude, as 
was done first in Refs. la. on smaller lattices, and 
more recently in Refs. [a, |9| on larger ones. The ensu- 



ing picture, which turns out to be closely related to the 
dual superconductor mechanism of confinement |lOj, es- 
sentially vindicates the scenario put forward by Einhorn 
and Savit 11] , in which the transition from the Higgs 
to the confined phase is triggered by proliferating vor- 
tices. As discussed below, the nature of this mechanism 
is consistent with the Fradkin- Shenker result 2] and the 
symmetry argument 4\ . Our results are at odds with the 
vortex string-breaking scenario put forward by Nagaosa 
and Lee ^3] j who recently argued that only the confined 
phase is present in the compact theory. Furthermore, we 
found no support for a deconfinement phase transition of 
the Berezinsky-Kosterlitz-Thouless (BKT) type recently 
proposed in Ref. 0], as for sufficiently large lattices we 
observe no scaling behavior at all across the phase bound- 
ary where the BKT transition is supposed to arise. 

The compact U(l) Higgs model is specified by the Eu- 
clidean lattice action S = S g + S<f,, with the gauge part 



Sg^P [1- COS 0^0)] 



(1) 



Here, (3 is the inverse gauge coupling, the sum extends 
over all lattice sites x and lattice directions /i, and 9^ v {x) 
denotes the plaquette variable 9 fiu (x) — A. lx 9 l/ (x) — 
A v 9 )Jl (x), with the lattice derivative A„# M (x) = 9^(x + 
v) — 9fj,(x) and the compact link variable Ofj,(x) £ [— vr, 7r). 
The matter part of the action S is given by 

S$ = -K,'^2p(x)p(x + fi) cos [A^x) - qOftix)] 

+ ]r{p 2 (.x) + A[p 2 (*)-l] 2 }, (2) 

X 

where polar coordinates are chosen to represent the com- 
plex Higgs field <f>(x) = p(x)e 1>fi ^ x \ with ip(x) G [— tt, it), n 
is the hopping parameter, and A the Higgs self-coupling. 
The pure |</>| 4 theory with fluctuating amplitude, ob- 
tained by taking the limit /3 — ► oo, was recently investi- 
gated by means of Monte Carlo simulations in Ref. . 
We consider the system in three spacetime dimensions, 



2 



taking one of the dimensions to represent (Euclidean) 
time. 

The precise nature of the phase diagram is investigated 
numerically by studying several (gauge-invariant) observ- 
ables chosen such that the gauge and matter parts are 
probed separately. Following Ref. 0, where the Lon- 
don limit of the model was considered, the gauge part is 
studied by measuring the monopole density M, as defined 
in Ref. 16], and the Polyakov loop. Both observables dis- 
tinguish a confined from a deconfincd phase. The matter 
part of the model is studied by measuring the Higgs am- 
plitude squared p 2 = (l/L 3 )J2 x p 2 (x), where L is the 
linear size of the cubic lattice. This bulk operator dis- 
tinguishes the Higgs phase from a disordered one. In 
addition, the plaquette action JTJ (divided by 3L 3 ) and 
the so-called coslink observable 

c = -^J2 C0 ^ A M*)-<iW} (3) 

are monitored. Both Metropolis and heat-bath meth- 
ods were used to generate Monte Carlo updates. Since 
these local updates become inefficient in regions of first- 
order phase transitions, the multicanonical method 01 
and reweighting techniques |l8j were implemented to ac- 
cess these regions of phase space. The simulations were 
carried out at fixed f3 on cubic lattices varying in size 
from 6 3 to 32 3 , in extreme cases to 42 3 . Thcrmaliza- 
tion of the production runs typically took 4 x 10 4 sweeps 
of the lattice, while about 10 6 sweeps were used to col- 
lect data, with measurements taken after each sweep 
of the lattice. The maxima of the coslink susceptibil- 
ity Xc — L 3 {(C 2 ) — (C) 2 ), and histograms, rescaled to 
equal height have been used to determine the location 
of the phase boundary. We have chosen the coslink sus- 
ceptibility to trace out the phase diagram because its 
peaks are more pronounced than for the other observ- 
ables. We have checked that within the achieved accu- 
racy, the monopole susceptibility peaks at the same lo- 
cation as xc does. Statistical errors were estimated by 
means of jackknife binning. For a detailed description of 
the algorithms and their implementation, see Ref. [l9| . 

The phase diagram 00, Q, summarized in Fig-.^for 
fixed gauge coupling (3 = 1.1, consists of two separate 
phases: a confined and a Higgs phase. In the lower right 
part of the phase diagram, the average Higgs amplitude 
squared takes on a minimum value (see Fig. [2J. The 
monopole density is finite here and practically indepen- 
dent of k and A. Snapshots of monopole configurations 
(see bottom inset in Fig. ^| show that the monopoles 
are in the plasma phase. As (3 increases, the monopole 
density decreases. The monopoles become completely 
suppressed in the weak gauge coupling limit (3 — > oo, 
where the model reduces to the pure |</>| 4 theory. The 
average plaquette action takes on a value also practically 
independent of k and A. These observables signal that 
electric charges are confined. This confined phase persists 




Figure 1: Phase diagram of the model at (3 = 1.1 in the 
infinite- volume limit. Solid dots mark the location of a line of 
first-order phase transitions. This line ends at a critical point 
around 0.031 < A c < 0.032. Open dots for A > A c mark the 
location of the Kertesz line (see text). Statistical error bars 
are smaller than the symbol size in the figure. The insets show 
snapshots of typical monopole configurations in both phases, 
with black dots denoting monopoles and grey dots denoting 
antimonopoles. 
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Figure 2: Averages of the Higgs amplitude squared p 2 , coslink 
observable C, plaquette action S g , and monopole density M 
as a function of k for various values of the self-coupling A, 
including the London limit A^oo, at/9 = l.lona relatively 
small lattice (L — 12). 



in the limit k — > 0, where the model reduces to the pure 
compact U(l) gauge theory first studied by Polyakov pj. 
Notice from Fig. |2 the non-monotonic behavior of (p 2 ) 
and (C) as a function of the Higgs self-coupling A. 

In the upper left part of the phase diagram, the aver- 
age (p 2 ) increases more or less linearly with increasing k 
(note the logarithmic scale used for p 2 in Fig. • The 
monopole density is vanishing small here and the few 
monopoles still present are tightly bound in monopole- 
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Figure 3: Maxima in the susceptibilities of the plaquette 
action S g , Higgs amplitude squared p 2 , monopole density M, 
and coslink variable C as a function of the lattice size L at 
P = 1.1 and A = 0.2. To fit all the data in one figure, Ixc** 
rather than Xc? ax i s plotted. 

antimonopole pairs 0] (see top inset in Fig. QJ. Being 
rendered ineffective, the monopoles can no longer confine 
electric charges. Taken together, these observables iden- 
tify this phase as Higgs phase. The identification of the 
two phases agrees with the behavior of the Polyakov loop 
we observed. 

We next examine how for given inverse gauge coupling 
(3 the confined phase goes over into the Higgs phase. Be- 
low a critical point A c (/3), we observe metastable behav- 
ior typical for first-order phase transitions, in accord with 
earlier Monte Carlo results obtained on smaller lattices 
0,0- Simulations for different values of f3 £ [1.1,2.0] 
show that the first-order phase transition becomes more 
pronounced for strong gauge coupling, i.e., small (3. At 
fixed [3, the transition becomes stronger with decreas- 
ing A, where fluctuations in the Higgs amplitude become 
more volatile For each value of A considered, the 
model was simulated on lattices of different sizes to study 
finite-size effects and to obtain precise estimates of the lo- 
cation of the first-order phase transitions, using the mul- 
ticanonical approach and reweighting techniques. The 
first-order line ends at a critical point, which for f3 = 1.1 
and infinite volume we estimated to be located in the in- 
terval 0.030 < A c < 0.032. We have not attempted to 
establish the nature of the critical point as critical slow- 
ing down requires very long runs of our code based on 
locale updates. 

Above A c , we observe the same remarkable behavior 
as previously found in the London limit A — > oo, where 
fluctuations in the Higgs amplitude are completely frozen 
0,0]. Namely, for sufficiently large lattices, the maxima 
of the susceptibilities do not show any finite-size scaling 
(see Fig. and the susceptibility data for the observ- 
ables in Fig. 2 obtained on different lattice sizes collapse 
onto single curves without rescaling, indicating that the 
infinite-volume limit is reached. Since a first-order phase 



transition can be excluded in this region, the absence 
of finite-size scaling suggests the absence of thermody- 
namic singularities in the infinite- volume limit. To also 
exclude a crossover in the usual sense, we have checked 
that the maxima of \c do not depend on the direction in 
which the phase boundary is crossed, either by varying 
A, or k |20j. This analysis is performed using multihis- 
togram reweighting techniques to achieve a high accuracy 
in determining the peak locations. To sum up this part, 
despite the presence of a well-defined and precisely lo- 
cated phase boundary, no ordinary phase transition or 
crossover in the usual sense seem to come into question 
above A c . 

In the dual superconductor scenario of confinement 
[lfT | , monopoles are pictured in the Higgs phase as being 
tightly bound together in monopole-antimonopole pairs, 
just as we observed (see top inset in Fig. QJ. The mag- 
netic flux emanating from a monopole is squeezed into a 
short flux tube, or vortex, carrying one unit 2ir/q (q = 1) 
of magnetic flux, which terminates at an anti-monopole. 
The vortices, which in this phase can also exist as small 
fluctuating loops, have a finite line tension. Upon ap- 
proaching the phase boundary, the vortex line tension 
vanishes. At this point, the vortices proliferate, gaining 
coiffigurational entropy without energy cost, and an in- 
finite vortex network appears which disorders the Higgs 
ground state. At the same time, the monopoles are no 
longer bound in tight monopole-antimonopole pairs but 
form, as seen in the bottom inset in Fig. ^ a plasma 
which exhibits charge confinement. This scenario ex- 
plains the existence of a well-defined phase boundary sep- 
arating the Higgs and confined phases [Jjj . The prolifer- 
ation of vortices persists in the weak gauge coupling limit 
f3 — > oo, corresponding to the pure |</>| 4 theory with global 
U(l) symmetry. Only in this limit, the proliferating vor- 
tices cause a continuous phase transition, belonging to 
the XY universality class in which the spontaneously bro- 
ken global U(l) symmetry is restored. Outside this limit, 
in the absence of a relevant global symmetry, the prolif- 
erating vortices do not lead to singularities in thermody- 
namic quantities and are not connected to a symmetry 
breaking transition. 

The situation is reminiscent of the Ising model in an ex- 
ternal magnetic field. It was shown by Fortuin and Kaste- 
leyn (FK) [2lJ that for zero field this spin model and its 
thermal critical behavior can be equivalently formulated 
as a correlated percolation problem by putting bonds be- 
tween nearest neighbor spins in the same spin state. The 
bonds are set with a temperature-dependent probabil- 
ity Pfk(T). The clusters thus constructed percolate pre- 
cisely at the Curie point and have the Ising critical expo- 
nents encoded in their fractal structure. By applying an 
external field, one explicitly breaks the global Z2 symme- 
try of the Ising model, and the partition function becomes 
analytic in temperature, excluding a thermal phase tran- 
sition. Yet, for a given applied field H, the FK clusters 
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still percolate at a precisely defined temperature T p (H). 
The resulting percolation line in the phase diagram is 
known as the Kertesz line 22]. Although percolation 
observables remain singular along the line, no thermo- 
dynamic singularities are encountered when crossing it 
[221 123|. In the limit H — > oo, all the spins are aligned 
along the field, so that the FK construction reduces to 
random bond percolation. The Kertesz line therefore 
ends at the temperature determined by pfk(T p ) = p c , 
with p c denoting the random bond percolation threshold. 
Along the entire Kertesz line, the percolation observables 
have the usual percolation exponents. 

The vortex proliferation line in the compact U(l) Higgs 
model is the analog of the Kertesz line in the Ising model 
in an external field. Such an interpretation of a decon- 
finement transition as a Kertesz line was firs t propo sed in 
the context of the SU(2) Higgs model [13, I^E E3l - Th e 
similarity with the Ising model can be made more precise 
by considering the London limit A — > oo of the compact 
U(l) Higgs model. In this limit, the vortex proliferation 
line starts at the XY critical point k — kxy,/3 = oo- 
To identify the endpoint, we note that for k — > oo, the 
model reduces to a Z q gauge theory. In 3D, such a dis- 
crete gauge theory is dual to the g-state Potts model 
with global Z q symmetry, which undergoes a phase tran- 
sition at some critical value q of the inverse gauge cou- 
pling 28] . Since the limit q — ► 1 of the Potts model de- 
scribes random bond percolation, we conclude that the 
vortex proliferation line ends at a random bond percola- 
tion critical point at (3 = (3\, k = oo. The vortex network 
present in the vicinity of this critical point is expected to 
be similar to the one studied in the context of the Kib- 
ble mechanism for cosmic string formation |29j, which 
was shown to belong to the random percolation univer- 
sality class 30]. By continuity, we expect that, although 
not connected to thermodynamic singularities, the per- 
colation observables have random percolation exponents 
along the entire vortex proliferation line. In the limit 
[3 — > oo, these exponents are expected to cross over to 
the ones appropriate for the XY universality class. In a 
future study, we plan to investigate the vortex network 
using percolation observables to numerically verify these 
conjectures directly. Such a study is more complicated 
than the indirect study of the vortex network presented 
here, via monopoles, which involves only the gauge sector 
of the theory and has the advantage that the relevant ob- 
servables (monopole density and corresponding suscepti- 
bility) can, in contrast to vortex percolation observables 
|3lLl32j| . be defined unambiguously. 

In conclusion, our Monte Carlo data on the 3D com- 
pact U(l) lattice Higgs model show that the Higgs and 
confined phases are separated by a well-defined phase 
boundary due to proliferating vortices. For fixed gauge 
coupling, the phase boundary is a line of first-order phase 
transitions at small Higgs self-coupling, which ends at a 
critical point. The phase boundary then continues as a 



Kertesz line across which thermodynamic quantities and 
other local gauge-invariant observables are nonsingular. 
In the London limit, the Kertesz line defined by the pro- 
liferating vortices connects the XY and random bond per- 
colation critical points, which both form limiting cases of 
the compact U(l) Higgs model. 
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